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04 ; I. INTRODUCTION, SUMMARY AND OUTLOOK 

O ; 

Regular lattices are used in statistical mechanics and field theory to describe crystal geometry or to discretize flat 
space-time. When dealing with irregular structures that arise in the study of many physical systems (for instance 
polymer gels, structural glasses or fractals) as well as in the discretization of curved space-time, more general structures 
■ are needed to describe the underlying geometry; prototypes of such structures are graphs, discrete networks made of 
nodes and links. 

t-H , In particular we are interested in "generic" , or even "random" realizations of infinite graphs, subject only to very 
general physical requirements. Consider, as examples in condensed matter theory, the infinite cluster of percolation 
theory or the fractal graphs of growth aggregates; other examples in more abstract contexts such as euclidean quantum 
gravity and random matrix theory, are graphs modelizations of both space-time and "target space" . 
\ Two physical requirements common to these applications are: first, the coordination of nodes is bounded (the 
I . number of nearest neighbors of an atom, molecule or basic building block has a geometrical upper bound); secondly, 
y—i ' in the limit of large radius the surface must be negligible with respect to volume (this is necessary if the structure is 
, embeddable as a whole in a finite-dimensional euclidean space) . 

Our main goal is the identification of suitable algorithms that can produce graphs which do fulfill those requirements. 
Moreover, since these graphs should be as "generic" as possible, that is with a minimal amount of constraints, a third 
natural property we ask for is that of "statistical homogeneity". By this we mean that the probability to find any 
finite neighborhood should be the same around every node. 

In this paper we restrict our attention to trees, that is graphs without loops, for two main reasons. They are 
important on their own, since they completely characterize special, physically relevant cases: for example, the incipient 
infinite cluster at critical percolation is essentially a tree |Q, while in two-dimensional quantum gravity interacting 
with conformal matter, for certain values of the central charge the metric collapses to that of branched polymers [^). 
Moreover, we want to use trees as a starting point for building graph: our project is to generate graphs by adding 
links to random trees and study properties such as the connectivity and spectral dimension as the density of loops 
increase. 

The outline of this paper is as follows. In next section we recall some definitions of graph theory, define an infinite 
bounded graph as a graph which satisfies the first two previously stated requirements, and recall the correspondence 
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between trees and branching processes. Next, in section III, we give a constructive definition of random trees by 



describing the algorithm used to build them and discuss their statistical homogeneity. In section IV we write the 
recursion for the growth probabilities and find by a simple and original method their scaling form in the limit of large 
sizes, from which we rederive the well known result that their local connectivity (or intrinsic Hausdorff ) dimension d c 
is equal to 2. Using statistical homogeneity we then obtain the same result for the average connectivity dimension d c . 

Along this way the following conjecture appears very natural: the average connectivity dimension for random trees 
is an upper bound for all bounded trees; that is to say d c < 2 for all bounded trees with saturation only in the case 
of random trees. There are indeed many examples of bounded trees with local connectivity dimension greater than 
two, for example NTo trees || or spanning trees of D-dimensional lattices with D > 2, but in these cases there are 
macroscopic inhomogeneities that cause d c to be different from its local counterpart d c and always such that d c < 2. 
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We think indeed that randomness and (statistical) homogeneity maximize the average connectivity dimension. 

II. BASIC DEFINITIONS 

We consider only infinite connected graphs and for any such graph Q we call G the set of its nodes (or sites, 
vertices). By definition G is in one-to-one correspondence with the set of natural number and any specific choice of 
such a correspondence defines an indexing, or labelling for Q . Unless otherwise specified we shall consider the graph 
as unlabelled, in the sense that we restrict our attention to properties that are independent on the labelling. 

The coordination (or degree) of a site x S G (that is the number of its nearest neighbors) will be denoted by z(x). 
The spherical shell S(x, r) and the spherical ball B(x, r) around any given site x are defined as 

S(x, r) = {y e G : d(x, y) = r} 
B(x,r) = {yeG: d(x,y) < r} 

where d(x, y) is the chemical distance on the graph; we shall denote the cardinality of these sets with s(x, r) — \S(x, r) \ 
(the surface area) and v(x,r) — \B(x,r)\ (the volume) respectively. In particular we evidently have v(x, 0) = 1, 
v(x, 1) = s(x, 1) + 1 = z{x) + 1 and in general 

v(x, r + 1) = v(x, r) + s(x, r + 1) 

The quantities s(x,r) and v(x,r) are example of "good observables" , in the sense that they describe properties of the 
unlabelled graph. More general observables of this type are the subsurfaces s z (x, r) and subvolumes v z (x, r) counting 
nodes with fixed coordination z in the shell S(x 7 r) and the ball B(x,r). 

The average of a function F : G — > R around the site x is defined as the infinite radius limit of the average on the 
ball with center x, namely: 

(F)(x)= lim — r^— J2 F ^ W 
i — >oo v(x, r) z — ' 

v ' y£B(x,r) 

Whenever (F)(x) does not depend on x, we drop the relative specification and regard it as the proper definition of 
graph-average (F). Then the measure of a subset A C G is identified with the graph-average (xa) of the characteristic 
function \A of the subset. In particular a central role is played by the the measures of the subsets of nodes with 
coordination z, namely the fraction of nodes with coordination z, which we denote by f z ; more precisely, we should 
consider first the limit 

f z (x) = km — (2) 

r^co V\X,r) 

of the fractions within a given ball and then worry about the dependence on x. We shall henceforth restrict out 
attention to /—graphs, that is graphs such that f z (x) does not depend on x, or f z (x) = constant = f z . By definition 
we have the normalization 



1 (3) 



Loosely speaking, one could say that f z represents the probability that a node chosen "at random" has coordination 
z. 

A. Local and average connectivity dimensions 

The quantities v(x, r) introduced above are sometimes called local growth functions of the graph. Their asymptotic 
behavior for r — > oo gives the (local) connectivity dimension d c [Q of the graph (sometimes called also intrinsic 
Hausdorff dimension) 

dc = i im lQ g;(^) 

r— >oo logr 
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The average growth function (v(r)) is the average of v{x,r) over the graph and its asymptotic behavior is related 
to the average connectivity dimension d c of the graph, defined as 

- \og(v(x,r)) 
a c = iim 

r^oo log?" 

Note that, even if it can be shown that d c does not depend on node x, d c can be different from d c because the limit 
r — * oo may not commute with the limit implied by the definition of the average. These two parameters are known 
to coincide on many "regular" graph (for example on .D-dimensional lattices) but are different on graphs which are 
manifestly inhomogeneous (as, for example the comb-like ^ and NTd trees ||). 



B. Bounded graphs 



We define an infinite connected graph Q as bounded whenever, for any x € G, the coordination is bounded 

z{x) < z < +00 



and the growth is volume-dominated, namely 



Km ^ = (4) 

r^oo v{X, r) 



Actually, it is easy to realize that if the second condition holds for one given node then by the first condition it holds 
for its nearest neighbors and then, by recurrence, for any other site. Moreover, these conditions imply that the average 
of any bounded function does not depend on the choice of the center of the ball [||. In particular, this means that a 
bounded graph is a /—graph. 



C. Trees and branching processes 



We shall from now on restrict our attention to trees, that is graphs without any loop. For a given connected tree 
one can follow a recursive procedure to find shells of increasing radius around any given site x: 

• the first shell is formed by all nearest neighbors of x and therefore its size coincides with the coordination of x: 

s(x, I) = z(x) (5) 

• for every node y of coordination z(y) in S(x, r) there is one link connecting y to a node belonging to S(x, r — 1) 
and z(y) — 1 connecting it to nodes in S(x,r + 1); since there are no loops, all links to S(x,r + 1) are directed 
to distinct nodes, so that 

s(x,r+l)= Yl My)' 1 } ( fi ) 

y£S(x,r) 

Thus the tree defines a specific realization of a branching process rooted at x, such that x has z(x) branches rooted 
on its nearest neighbors while every other node y on S(x, r) has z(y) — 1 branches rooted on its nearest neighbors on 
S(x, r + 1) (notice that we may equivalently consider all branches as rooted on links rather than nodes). In turn this 
identifies a coordination sequence {z\, Z2, 23, . . . }, where z\ = z(x) is the coordination of x, Z2, 23, . . . , z v r Xt i) are the 
coordinations of nodes on the first shell, zw^n+i, z v (x,i)+2i ■ ■ ■ > z v(x,2) are the coordinations of nodes on the second 
shell, and so on. The indexing on each shell is uniquely fixed by the branch structure as soon as the branch roots on 
each nodes are indexed (we assume in a consecutive manner). Therefore a coordination sequence identifies a unique 
labelled and rooted tree, while there are many coordination sequences, obtained by permuting labelled branches over 
any given node and by changing the root itself, that correspond to the same unlabelled tree; the latter identifies 
therefore an entire equivalence class of coordination sequences, or isomorphic labelled trees. 
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D. Bounded trees 

Given a branching process that builds an infinite rooted tree, from equation (^) one can obtain the rate of growth 
of the shells 

s(x,r) 

s(x, r + 1) - s(x, r) = ^2 [z{y) - 2] = ^ (z v(x . r - 1)+j - 2) 

»eS(i,r) 2=1 

which, after integration and use of equation (||), gives 

v(x,r) 

s(x,r + l) = 2+ ]T [<»)"2]= 2+ £ (^-2) 

yeB(x.r) j=l 

In the limit of infinite radius, if the tree is bounded according to equation (|j) , the average coordination follows 

, , six, r + 1) , . s(x, r) 

(z ) - 2 = lim y I ' <(z-l) lim ) - = 

We have dropped any reference to a; in (z) since this quantity does not depend on it for a bounded graph. For this 
reason, by considering shells centered around any other node x 1 of the tree rooted at x, the same average coordination 
(z) = 2 would have been obtained. Thus (z) = 2 is a property of the unrooted tree. In fact, since a bounded tree is 
an /—tree, it may equivalently be written as an average over the "probabilities" f z : 

(*> = X>/, = 2 (?) 



E. Link orientation 

We have seen that a rooted tree has a natural orientation for all links determined by the branching direction: the 
root x has z(x) outgoing links and any other site y has z(y) — 1 outgoing links and one incoming links. Thus, with 
the unique exception of the root, there is a one-to-one correspondence between nodes and incoming links, so that 
f z is also the fraction of links leading to nodes with coordination z; on the other hand the correspondence between 
nodes and outgoing links is one to z — 1 so that the fractions of links departing from nodes with coordination z is 

f z = (z- l)f z (8) 

Note that the f z may be interpreted as properly normalized probabilities on a bounded tree because of equations (|^) 
and ©: 

Z Z 

E/. = E(*-i)/* = 2-i = i 

2 = 1 2 = 1 

III. RANDOM TREES 

There exist many different ways of defining and generating random trees (random binary trees, pyramids of various 
order, linear recursive trees, branched polymers and many others). We restrict our attention to growth stochastic 
algorithms that are local and produce infinite trees that are bounded and statistically homogeneous (see below). 

A. Random branching processes 

Taking into account the correspondence between trees and equivalence classes of coordination sequences {zi, Z2, ■ ■ ■ } 
as described above, the most immediate definition of such a random tree with predefined coordination fractions 
{fz, z = 1, 2, . . . , z} fulfilling equations (0) and (0), would be as a random coordination sequence in which each Zj 
is extracted from the set {1,2,... ,z} with probability f Zj . This is known as the (critical) Galton- Watson stochastic 
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branching process : it satisfies the locality requirement that the coordination of one node does not depend on that 
of any other node and guarantees that the average coordination is two, but does not fulfill our request that the tree 
grows indefinitely (the growth stops when the last shell of the tree is made solely of node with coordination 1). In fact 
a classical problem of probability theory is the estimate of the survival probability of such a branching process: one 
finds that the probability that a tree has more than N nodes to be of order iV -1 / 2 0. Another shortcoming of this 
approach, as presented above, is the bias on the tree root due to fact that it has as many branches as its coordination, 
while every other node has one branch less than its coordination. This is easily remedied by constraining the tree 
root to have coordination one; one must then consider coordination sequences of the form {zi, Z2, 23, ... } with z\ = 1. 
Then all nodes with index greater than 1 are statistically equivalent in the sense that: i) each one of them is the root 
of a new random branching process distributed exactly as the one rooted on node 1; ii) processes rooted on different 
branches are statistically independent. 

To overcame the extinction problem, one may consider an infinite coordination sequence as a collection of dis- 
connected finite trees of any size (sometimes called random forests or branched polymers ), by simply using the 
first coordination extracted after extinction as first branch root of new branching process. This corresponds to a 
"grand-canonical" ensemble for the statistic of random trees [§. 

We are looking instead for a process preconditioned on non-extinction @] that builds one infinite connected bounded 
tree "at random" . Of course, for practical numerical purposes, to build very large connected trees one needs only 
to make several attempts with the standard Galton- Watson branching process (with the root coordination z(x) con- 
strained to be one or extracted with probability proportional to f z /z) until the required size is reached. Alternatively, 
to ensure that the tree grows indefinitely we must implement in the branching process itself the condition of non- 
extinction (this is our choice) . The link orientability play here a crucial role. Suppose in fact that an infinite connected 
tree is already given and that we reconstruct it starting from a certain node x by extracting the right coordination 
sequence; x is naturally the root of the corresponding branching process, but since we assumed the infinite tree to be 
there already, we may choose another node o as root from which the link orientation propagates; for arbitrarily large 
reconstructed balls B(x,r), we may assume that o does not belong to B(x,r), that is r < d(x, 6) (there is a set of 
measure one of candidates o which satisfy this condition on an infinite tree); consider now, for any r' — 1,2,... , r, 
the s(x,r') links between S{x,r') and S(x,r' — 1); one and only one of them is incoming, that is oriented towards 
the inner shell, while the other s(x, r') — 1 are outgoing, that is oriented towards the outer shell; taking equation 
into account this is the essential observation to define the stochastic building algorithm we are looking for; we call it 
algorithm A: 

Al: extract the coordination z\ for the root node x with probability f Zl , then pick at random, with equal probability 
1/zi, which of the Z\ link on x is the incoming one and mark it; 

A2: on the shell S(x, r), for r = 1,2,3,..., extract the coordination zm of the node attached to the marked link 
with probability f ZM , then pick at random, with equal probability (zm — 1) , one of the newly added links and 
move the mark on it; extract the coordinations Zj of the other s(x,r) — 1 nodes on the shell with probability 

JZj ■ 

We see that the marked links are exactly those forming the unique path connecting x to o. As this stochastic 
algorithm proceeds indefinitely, o is removed to infinity. Since the random choosing events in the algorithm are all 
independent, the probability of a coordination sequence corresponding to B(x,r) reads therefore (here v — v{x,r)) 



Prob[-{>i,32,... z v }\ = — — TT f z{y) 

s(x, 1) J - J - ' 



s(x, r + 1) 

11 Iz{y) 

(9) 



1 



E 

3=1 



where, except for the first shell, all factors of the form (zm — l)^ 1 have canceled against the factor zm — 1 hi fz M (see 
equation dq)), and the factor s(x,r + 1) comes from the sum over all possible positions of the marked link on the last 
shell S(x, r) (this is indeed the probability for an unmarked tree). We see that non-extinction is realized by forcing at 
least one node per shell, the one which has the marked link as incoming link, to have a coordination z > 2 extracted 
with the modified probability (^). This is done however in such a way that no trace of the marks remains except that 
for the last completed shell S(x, r), which is eventually removed to infinity; for finite B(x, r) even this residual mark 
disappears by summing over all its possible positions, as in expression (^|). 
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B. Statistical homogeneity 

The probability (|^) is clearly invariant under all permutations of the coordinations that do not touch Z\. The 
different role of the coordination Z\ — s(x, 1) with respect to all the others is due to the fact that coordination 
sequences correspond to rooted and labelled trees. This asymmetry disappears when we consider the tree as unlabelled 
and therefore multiply (||) by the number of distinct labellings of the rooted B(x,r), This number is proportional, 
through a suitable symmetry factor, to Z]\ IX^i^j — -0' smce the root has Z\ branches, while all other nodes have 
Zj — 1 branches. By the same token, we could consider other algorithms for producing coordination sequences as long 
as they provide the same probability for the unlabelled B(x,r). It is not difficult to realize that one such algorithm 
is the following algorithm B: 

Bl: extract the coordination z\ for the root node x with probability / Zl ; 

B2: on the shell S(x, r), for r — 1,2,3,..., extract the first coordination z — 2W X)r _i)+i with probability / z , then 
extract the coordinations Zj of the other s(x,r) — 1 nodes on the shell with probability f Zj . 

With this algorithm a specific path is explicitly selected on the tree, but this leaves no traces on the infinite random 
unlabelled tree. In more precise probabilistic terms the situation is described as follows: 

• every instance of an infinite critical Galton- Watson branching process conditioned on non-extinction has almost 
surely a unique path T (the spine) extending from the root to infinity Q . 

Evidently the spine is just our path of marked links. As a consequence, there exists the so-called spinal decomposition 
of the infinite tree |9| , which in our case, with the root coordination z\ extracted with probability f Zl , may be stated 
as follows: 

• the branching process that builds the tree may be regarded as a collection of infinite independent subprocesses; 
one processes builds the spine T by extracting the coordinations of its nodes with probability f z ; each node 
y G T is the root of z(y) — 2 (z(x) — 1 if y = x is the root) independent and identically distributed critical 
Galton- Watson processes without condition on non-extinction; 

We see therefore that if the tree is cut on every link of the spine T, a Galton-Watson random forest is obtained. As 
discussed above, this "grand-canonical" ensemble of random trees is statistically homogeneous. It is also clear that 
the spine does not induce any statistical inhomogeneity on the infinite connected random tree, since all traces of the 
marks which identify the spine disappears as we have already noticed in the calculation of the ball probability (see 
previous section); moreover, as we shall see in the following, the spine is a zero measure set in the tree. 

IV. GROWTH ON RANDOM TREES 

In this section we study the properties of local and average growth functions on the infinite random trees introduced 
in the previous section. Let us begin with the growth around the node x used as root in the building stochastic 
algorithm. 

A. Recursion rules 

From the explicit form (^|) of the probability we may quite easily derive recursion rules for the probability distribu- 
tions of the random variables v(x, r) and s(x, r) (volume and surface area of the ball), that is 

P r (w, s) = Prob [ v(x, r) = v and s(x, r) = s ] 

To denote the probability distribution for v or s only, we define 

S 

P r (>,8)=Y,Pr(v,s) 

v 

When r = 1, from equation (||) we read the probability of a ball made of the root only, so that 

P 1 (v,s) = f s 5(v-s-l) 
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where S(n) = 1 if n = and zero otherwise, while the use of the iterative rule gives the Markovian recursion 

Pr+lW, S') = S ) 5 ( V> -J'") (10) 

v,s 

with the surface-to-surface transition probabilities, easily derived from (^|), 

^-' = 7 E mX> -!)-*') II/** (ii) 

21,2:2,... ,z s \i=l / i— 1 

Now it is useful to introduce the generating functions, or discrete Laplace transforms 



oo oo 



G r (A^)=^]TP r («, S )AV 



V — 1 s— 1 



for which the recursion reads 



G r+ i(A, M ) = E E W«V)aV 

u' = l S ; — 1 

OO OO oo 

= E E p >> s ) A " E (W 

u— 1 s — 1 s' — 1 

oo oo 

= E E P '( v > s)X"X^g'(X fJ ,) giX^- 1 



v—l s—1 



that is 



where 



G r+1 (A, /1 ) = A^G r (A,j(A fl )) (12) 



5 (A) = (A 2 - 1 ) = J2 = h + / 2 A + / 3 A 2 + . . . 

2=1 

is a polynomial in A with the properties 

g(l) = l , g'(l) = l (13) 
which follow from equations (||) and (0). At r = 1 we have 

Gi(A,/i) = \ 2 ng(\n) 

while for A = /i = 1 we verify 

G r+1 (l,l) = G r (l,l) = G 1 (l,l) = l 

as required by the proper normalization for P r (v,s). By construction G r (l,/i) is the discrete Laplace transform of 
P r (- , s) and satisfies the simpler recursion 

G r+1 (l )M ) = /i^G r (l, 9 (/i)) (14) 
9W 

which may also be written as 

G r+1 (l,n) = fj,g r+1 (n)g' r (n) , g r+ i(n) = g(g r (n)) , So(m) = M (15) 
Similarly G r (A, 1) is the discrete Laplace transform of P r {v, •) and satisfies the recursion 

G r +i(A,l) , G r (A,l) 

K +1 (x) =X9{hr - lix)) -hAxT (16) 
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where h r (X) fulfills its own independent recursion rule 

h r+1 (\) = \g(h r (X)) (17) 

with the initial condition /iq(A) = A. 

It should be stressed that the recursion rules ©, (0) and @ fr om a numerical point of view represent already a 
complete solution of the original problem of determining the probabilities P r (v, s), P r {- , s) or P r (v, •) for any r (when 
r is small the problem is almost trivial for any symbolic manipulation package). In fact we can choose A and /x to be 
the roots of unity of order Ni and -/V2, respectively, with Ni < ?) max and N2 ^ s max , where f max , s max are the largest 
values at which P r (v,s) is nonzero. Then the recursion rules are used to numerically find G r (A,/i) from the initial 
condition and finally P r (v, s) is approximated by inverse discrete Fourier transform from G r (A, /i) (of course, for large 
r we must take N± and iV 2 large but fixed, much smaller than w max and s max ). We have performed this program 
separately for surface and volume probabilities using the numerical package Matlab with the results plotted in figures 
I and f| 

B. Expectation values and scaling 

The expectation values for the volume and surface area, defined as 



S 

E(v(x, r)) = vPr(v, ■) = ^G r (X, 1) 



fj.=l 



can be explicitly calculated from the recursion rules; in fact for the first shell one simply gets 

E(s(x, 1)) = (z) = 2 
E(u(a:,l)) = 3 
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FIG. 2: Numerical evaluation of P r {v, 



while for larger values of r the recurrence relations are easily found to read 

E(s(x, r + 1)) = E(s(x, r)) + 2a 

and, quite obviously 

E(v(x, r + 1)) = E(w(x, r)) + E(s(x, r)) + 2a 



where 



Thus we obtain 



2a = g"(l) = (z 2 )-4=(z 2 )-(zf = (Az) 2 



E(s(x, r))=2 + 2a(r - 1) , E(v(x, r)) = 1 + 2r + a(r - l)r 



(18) 

By further differentiating equation (|l2|), one finds recursion rules for the expectation values of higher powers of s(x, r) 
and v(x,r) as well as of product of them (the moments of P r (v, s)). By differentiating the logarithm of ( p^|) one gets 
the centered, or connected moments. For instance, the general form of the connected surface moments is found to be 



E c ([s(a;,r)] fc ) = c k a k r k + 0(r k ~ v ) 



(19) 



with c/c constants; this fact suggests the scaling hypothesis that in the limit r — > oo of the shell probability must be a 
function only of the ratio s/(ar) (this is confirmed by the numerical determination in figure 0, after proper rescalings) 



arP r (- ,s) 

where the factor ar on the left comes from the integration measure in the normalization: 

1 = J2-(arP r (-,s)) r ~°° / dr0(r) 



(20) 
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By the same token we see that the functions 

CO ^ 

F r (u) = G r (l,exp(-^)) = Yl — (arP r (-,s)) exp(- 
should converge as r — > oo to the (continuous) Laplace transform -F'(u) of the the scaling function </>(t) 

/>00 

F r (u) r -°° / dr0(r)e- T " = F(u) 



To calculate we start from equation (15) and notice that the property ( |13| ) and the recursion g r+ i(/j,) = g{g r {^)) 
imply g r (l) = g' r (l) = 1 for all r. This entails the scaling law 

ffr ( exp (-^)) =i + ^M + fi 

with A(0) = and A'(0) = -1. Therefore 

F r (u) = -arg r (exp(-^))-^-g r ^(exp(-^)) ~-A'(u)+o(l) (21) 

yielding the identification F(u) — —A'(u). Then the recursion rule for g r (p) to leading non-trivial order becomes the 
differential equation for A{u) 



with the solution 



Hence 



uA' = A(l + A) 
1 + u 



(1 + u) 



with the inverse Laplace transform 

so that the surface probability reads for large r 



(r) =Te~ T 



P r (-,s)~ —e~ s l^ (23) 
a r 

We may repeat this approach for the volume probabilities as well, by assuming a scaling form in the variable v/(ar 2 ) 
(this would be obvious from the surface scaling for an homogeneous structure and is confirmed by the first connected 
moments and by the numerical data of figure |2j). Thus we set (notice that g(l) — 1 implies h r (l) — 1, see equation 
©) 



Mexp(-^))=l + ^ +0 (i) 



Substituting this expression in equation ([l7]) we obtain to leading non-trivial order 

2£B' = B(l + B) - £ 

whose solution, taking care of the condition B'(0) = —1 that comes from h' r (l) — r, reads 

B(0 = -Vetanh(70 

Now we set 

G r (exp(-^),l)=ff(0+o(l) 
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so that 



G r+1 (exp(-^),l) ^ 2^^ Jl\ 
G r (exp(-^),1) " + r + 

and since for A = exp(— jjs) 

A5W)) w" 1+ ? B(6+o 

substituting these expressions in equation ( |l6| ) we obtain the equation for 
and finally 

ff <« = < 24 » 

due to the normalization condition H(Q) = 1. The inverse Laplace transform of -ff (£) is just the scaling form of the 
volume probabilities 

and may be written as (see Appendix [A] for details) 

$(r) = m 1 ' 4 f^^V 2 (^(1 - m)^(m) - l) 

where 

K(l - to) 



and K(m), E(m) are the two complete elliptic integrals 



/•t/2 J/) /■t/2 

K{m)= I —^^^=, E(m)= d6Vl-msm 2 t 
Jo v 1 — msin 2 f Jo 



The large r behavior of $(t) follows in the limit to — > 0, 

to - 16 

and reads 

$(r) ~ (tt 2 t - 2) e -*V4 , T -> oo 
The behavior for small r is related to that for m — ► 1, 

1 - to ~ 16e~ 1/T 

so that 

$( T )~ JL T -3/2 e -l/T r ^ 

Therefore the volume probability distribution has the large r asymptotic behavior 

„-ar 2 /v for „ < r 2 



(25) 



7rl/2 w 3/2 

Pr(v,-)^{ 2 (26) 
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C. Local and average connectivity dimensions 



The results of the previous section concern the probability distributions for the growth around x in the statistical 
ensemble of infinite trees rooted at x. To relate these results to graph averages over a single "generic" realization of 
such a tree, we consider the following observables 



Pr(v) 



1 



v(x, R - 



r) ^ 



S(v(y,r) - v) 



that are the fractions of nodes y for which the ball of radius r has volume v in a specific realization of a tree rooted 
at x (restricted to radius R). In the limit of R — > oo the natural ergodic, or better autoaveraging hypothesis is that 
P r (v) deviates from its expectation value E(P r (u)j only on a subset of measure zero of the statistical ensemble and 
that E(P r (w)j coincides with P r (v, •), the previously calculated probability 



lim AP r (v) = 

R — ►qq 



lim E(P r (u)) = P r (v,-) 



(27) 



As a matter of fact, the results reported in section III B are enough to conclude that this autoaveraging is not simply 



an hypothesis: let us fix the size r of balls B(y, r) centered on arbitrary nodes and consider the local finite branching 
process rooted on y that reconstructs P(y,r); if y dists more than r from from the spine this process is distributed 
exactly like that rooted on x; if two such processes are rooted on nodes well separated, they are independent; the 
number of nodes closer than r to the spine scales with i?, the size of the main process, as compare to the volume 
v(x, R — r) that scales as R 2 ; from the (weak) law of large numbers we then expect ( p7) to h old. We have performed 
several numerical checks that confirm this expectation with high accuracy (see section IV D| ) 

We conclude therefor e, b y looking at the main results (p6|), that the local connectivity (or intrinsic Hausdorff) 
dimension d c of section [I A exists and that d c = 2 with probability one. Moreover, d c = d c because, thanks to (|27 
we can substitute the expectation value E(v(x, r)) in Jlg| ) with the graph average (v(r)) for a single random tree. 



D. Numerical checks 

We have performed several numerical checks of the autoaveraging property discussed above; for example figure |3| 
compares the distribution of r = v(y,r)/ (ar 2 ) for 3 • 10 5 nodes y on a graph of 10 7 nodes with r — 500, with the 
calculated limit distribution $(r). 

Another check can be done in the case of v z (x, r), since its average over different realizations of the graph-generating 
algorithm can easily be calculated and compared with the numerically obtained average over different nodes on a single 
graph. The spinal decomposition and the considerations in the previous section allow us to write 

E(«,(x,r)) = /,r+(E(t;(s,r))-r)/, 
= f z (l + (z - a)r + ar 2 ) 

with the use of equation ([l8|). In figure ^ we have plotted E(w z (x, r)) / (fzE(v(x, r))) compared with the corresponding 
average over 3 • 10 4 nodes of a single random tree with 10 7 nodes. 



APPENDIX A: INVERSE LAPLACE TRANSFORM OF 



The inverse Laplace transform of H(£), $(r) is defined to be 

p7+i oo 



$(r) = — / d^(0exp(r0 



where 7 is an arbitrary positive constant chosen so that the contour of integration lies to the right of all singularities 
in H(£). In this case H(£) is given by equation (M) 
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0.5 1 1.5 2 2.5 3 3.5 4 4.5 



T = v/(oc r 2 ) 

FIG. 3: Comparison between the distribution of v(y,r)/r 2 (r — 500) in a graph and $(r) 




FIG. 4: Comparison between numerical averages over 3 • 10 4 nodes of a single tree with 10 7 and expectation values. Error 
bars are smaller than the size of symbols 
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and its singularities in the complex plane of the variable £ are located in 

C„--7r 2 (n+i) 2 n = 0,1,2,... 
and the residue of H (£) exp(r£) in £„ reads 

Res[if (£) exp(rO, f„] = (tt 2 (2t7 + l) 2 r - 2) e -" 2 («+i/2) 2 r 
Now using the residue theorem 

oo 

*(r)=£Rea[fr(Oexp(rO,Cn] 



= J2 (^i 2 ™ + l f T - 2 ) e^^ n+1 '^ T ( A1 ) 
n=0 

= (-2r|--l)^(0,e-^) 

where 62(11,, q) is the second elliptic theta function with argument u and nome q; for it = it can be expressed in 
terms of the complete elliptic integral K[m) 

with m obtained by inverting 



Using the identities 



and 



we obtain 



K(l - m) 

T = 

nK(m) 



d E(m) - (1 - m)K(m) 

dm {m) ~ 2m(l-m) 



K(m)E(l - m) + K(l - m)E(m) - K(m)K(l -m) = - 



8 d 
-It— — —to(1 — m)K(m)K(\ — m)— — 
or 7r dm 



and finally equation 



$(r) = m 1 / 4 (^^) V2 ~ m)^(m) - l) 

Now the asymptotic behavior for r — > oo and r — > can be found using the expansions 

#(m)~-(H ) for 

. 7T / 777, N 

E.(m) ~— (1 — — J tor 777 — > 

K(m) ~ log 4 — — log(l — m) + — (1 — m) (log 4 — — log(l — m) — 1 ) for m — > 1 



J5(m) ~1 + -(1 - 777 ) (2 log 4 - log(l - 777) - 1) for 777 — > 1 



For to — > one can obtain 

4 log 2 — log m 

T ~ ^ 

to ~ 16e _7r2T 

so that the asymptotic behavior $(r) for r — > +oo 

$(r) ~ to 1 / 4 ( log 4 - i log m - l) 



tar that is simply the first term in the series equation ( ]Al| ) . 
For to — > 1 the variable change reads 



1 



" 4 log 2 - log(l - to) 
l-TO~16e~ 1/r 



so that for r — > 



*(r) ~ (4 log 2 - log(l - to)) 3/2 



4 T -3/2 e -l/r 
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